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O ■ Abstract 

^ ■ We discuss the Pade approximation to the tttt scattering ampli- 

q^. tudes in 1-loop chiral perturbation theory. The approximation re- 

D | stores unitarity and can reproduce the correct resonance poles, but 

' the approximation violates crossing symmetry and produce spurious 

CN ■ poles on the complex s plane and therefore plagues its predictions on 

| \ physical quantities at quantitative level. However we find that one 

' virtual state in the IJ=20 channel may have physical relevance. 
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! The chiral perturbation theory is a powerful tool in studying strong 

Q_i! interaction physics at low energies. The principle of chiral perturbation the- 

p ' ory is that it incorporates the global symmetry of the QCD Lagrangian - 

the spontaneously broken chiral symmetry - into the effective Lagrangian, 
and directly deals with the physical states - the pseudo-Goldstone bosons 
- as the physical degrees of freedom in the effective Lagrangian. The phys- 
ical amplitude is calculated by a perturbative expansion in powers of the 
light quark masses and the external momentum, and unitarity is respected 
perturbatively. As the external momentum increases, however, the chiral 
expansion diverges rapidly and the unitarity condition is badly violated. To 
remedy such a situation, the Pade approximation has been used to improve 
the behavior of chiral amplitudes at higher energies, which has stimulated 
revived interests in the recent literature (see for example 0-0). One of the 
cost in using Pade approximation is the violation of crossing symmetry which 
has also been discussed in the recent literature (see for example Refs. [0, || 
andi). 
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For the ttty — > 7T7T scattering process, the isospin amplitudes in the s 



channel can be decomposed as [|10 

rpl = 



T 



1=1, 
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T I=2 (s,t,u) 



3A(s, t, u) + A(t, u, s) + A(u, s, t) 
A(t, u, s) — A(u, s, t) , 
A(t, u, s) + A(u, s, t) , 



(1) 



as a result of the generalized Bose statistics and crossing symmetry. In 
SU (2) x SU(2) chiral perturbation theory (ChPT) to one loop || we have, 
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+15m* + 12m 2 (s - m 2 )^ - 3m%} 
93.3MeV is the pion decay constant and the function J(s) is 
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In Ref. the U parameters are taken to be 
h = -0.62 ± 0.94 , l 2 = 6.28 ± 0.48 , l 3 = 2.9 ± 2.4 , l 4 = 4.3 ± 0.9 (4) 



determined from low energy experiments. 

The partial wave expansion of the isospin amplitudes is written as 

T\s,t,u) = 32tt^(2J+ l)Pj(cosfl)Tj(s) , 
j 



(5) 



where, for titt scatterings, the sum is over even(odd) values of J for even(odd) 
values of / because of the restriction of Bose statistics. The inverse expression 
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047T J-i 

2t 

cosO = 1 + -— , (6) 

s - Ami 

u = Ami — s — t . 

The partial wave amplitudes in ChPT expanded to 0(p 4 ) are, 

T I As)=Ti 2 {s)+T I JA {s), (7) 

where Tj 2 (s) and Tj 4 (s) represent terms of order 0{p 2 ) and 0(p 4 ), respec- 
tively. The Tj 2 (s) amplitudes can be rigorously derived from current algebra 
and are model independent. The [1,1] Pade approximation to the partial wave 
amplitudes is, 

From perturbative unitarity in ChPT we have 

ImTj j4 ( S ) = p\T^ 2 (s)\ 2 , (Ami < s < 16ml) . (9) 

With this relation, it is easy to prove that the [1,1] Pade approximants given 
in Eq. (§) satisfy elastic unitarity: 

lmT I /' 1 \s)=p\T I J ^(s)\ 2 , (10) 

and it is well known that the phase shifts they predict are considerably 
improved comparing with the perturbation results. However, the total am- 
plitude 

T(s, t, u) = 327r^(2J + 1)Pj(cos0)Tj{s) , (11) 
j 

has certain crossing properties which are lost in the amplitude defined by 

T p (s,t,u) = 32vr^(2J + 1)Pj(cos6)t}' 1] (s) , (12) 
j 

by violating the so called Balachandran-Nuyts-Roskies relations, as discussed 
recently in Refs. |7], ||. 
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The Pade approximation not only gives an improved prediction to the tttt 
scattering phases comparing to the perturbative results, but also it predicts 
the correct physical resonances on a qualitative level. For example in the 
amplitude Tq' 1 ' 1 ^), on the second sheet of the complex s plane one finds the 
u resonance with M a = 430MeV and IV = 456MeV, [] and the p resonance 
in the IJ=11 channel: M p = 708MeV and T p = 119MeV which should be 
compared with the experimental value M p = 769MeV and T p = 150MeV. 
Actually, the quality of the Pade approximants can be much improved by 
tuning the U parameters within reasonable ranges, that is to fit the global 
phase shifts while the U parameters are still consistent with the parameters 
constrained by low energy experiments. This remarkable improvement is 
possible, as pointed out in Refs. f|, H- Inspired by this approach (often 
called unitarized chiral approach in the literature), we also made a global fit 
to the experimental phase shifts using the [1,1] Pade amplitudes in Eq. (|8|). 
In our fit the experimental data are taken from Refs. [|T2J especially 



we also include the newest E865 Collaboration data from K ei decays Jl6| 
We fit the IJ=11 and 20 data up to \fs = lGeV and the IJ=00 data up to 
y/s = 700MeV. The fit quality is impressive, as can be seen in Fig [l|. The 
resulting values of the U parameters corresponding to the minimized \ 2 are 
as follows, 

fi = -0.485 , l 2 = 5.696, l 3 = 16.136 , Z 4 = 3.650 . (13) 

Notice that in the above fit the 1$ parameter is much larger than that in 
Eq. ([|) and the value in Ref. ||. This is not a serious problem as we found 
in the fit that though the minimization of the y 2 prefers a large value of Z3 
but it is not sensitive to I3. Another fact which is worth pointing out is that 
the chiral estimate on l 3 is crude and is subject to large uncertainties |T7 



The physical resonances obtained from our global fit are, M p = 751MeV, 
r p = 144MeV and M a = 449MeV, r CT = 482MeV. The scattering length 
parameters in the IJ=00 and 20 channels are found to be a(] = 0.194 and = 
—0.044 which are reasonable as comparing with the experimental results. 

Though successfully predicting the existence of the a and the p resonances 
qualitatively or even quantitatively for the latter, it is well known that the 

Using the same set of U parameters as in Eq. (g) and the method of Ref. 0, a fit to 
the experimental phase shift in the IJ=00 channel, with the left hand cut (l.h.c.) obtained 
from Eq. (§), gives instead M CT = 519MeV and = 579MeV. 

2 The data truncation in the IJ=00 channel is to reduce the coupled channel effects and 
the pollution from / (980). 
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Figure 1: The global fit in IJ=00, 11 and 20 channels, using the [1,1] Pade 
approximants and the procedure described in the text. 
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Table 1: Resonances, physical sheet resonances (PSR), bound states (BS) 
and virtual states (VS) as predicted by Pade approximation on the complex 
s plane using the U parameters from Ref. ||. The pole position s po i e = 
(M + iT/2) 2 . Here the residue of a second sheet pole means the residue of 
the corresponding pole in 1/S. 



Pade approximation encounters severe problems in theory aspects. One of 
which is the generation of dubious poles on the complex s plane, as listed in 
tables 1 and 2, | except for those physically accepted. We categorize these 
dubious poles into 3 classes: the first contains the bound state pole and the 
accompanying virtual state pole denoted as BSi and VSi respectively (the 
latter corresponds to the zero in the physical sheet below threshold) in the 
IJ=00,20 channels. The positions of the pole and the accompanying zero are 
very close to each other. The second includes the resonance poles found in 
the physical sheet (denoted as PSR), and the third class is for the virtual 
state (denoted as VS 2 ) near s = found in the IJ=20 channel. 

There are two reasons to ignore the first class poles. The first reason is 
that they have very tiny residues and can be safely ignored numerically. The 
second reason is that they can be tuned to disappear totally by varying the 
U parameters in the reasonable range. To see this more clearly we recall that 

3 The poles listed in tables 1 and 2 may still be incomplete. 
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Table 2: Resonances, physical sheet resonances (PSR), bound states (BS) 
and virtual states (VS) as predicted by Pade approximation on the complex 
s plane, using the k parameters obtained from the global fit to the phases 
in the IJ=00,20,11 channels as described in the text. The pole position 
s poie — (M + iT/2) 2 . Here the residue of a second sheet pole means the 
residue of the corresponding pole in 1/S. 



the S matrix of the [1,1] Pade approximant takes the following form,Q 

o[i,i] = T 2~ T 4* = T 2 -ReT 4 + zp(T 2 ) 2 

" t 2 _t 4 T 2 — ReT 4 — ip(T 2 ) 2 ' 1 ) 

where the spin and isospin indices are dropped for simplicity. A virtual 
state is an S matrix zero located on the real axis below the threshold which 
requires T 2 = T 4 whereas the bound state pole requires T 2 = T 4 . Therefore 
the bound state pole and the virtual state disappear simultaneously if the 
following requirements are met, 

ImT 4 = , ReT 4 = T 2 . (15) 

Since ImT 4 = p(T 2 ) 2 from perturbative unitarity the first condition in Eq. ( |15|) 
implies that the bound state pole and the virtual state pole cancel each other 
when they move towards the Adler zero position for the lowest order partial 
wave amplitude. As a consequence, the Adler zero of the Pade approximant 

4 In the following all terms should be understood as their proper analytic continuation 
on the complex s plane. For example, ImT 4 = p(T2) 2 and the r.h.s. is well defined on the 
complex s plane. 
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is still there which coincides with the one of the lowest order partial wave 
amplitude, but the order of the zero decreases from two to one (the sec- 
ond order Adler zero is unnecessary from theoretical point of view). The 
second condition of Eq. (|T5D, which now reads ReT 4 = 0, affords a con- 
straint that those /, parameters have to obey. f\ To cancel the bound state 
in the IJ=20 channel, we can for example choose the following set of If. 
Ti = —0.5,l 2 = 5.3, T 3 = 3.7, U = 4.0, satisfying Eq. (|i~6|) numerically which 
fit the experimental values of m p and V p well. For these U parameters the 
experimental phase shifts of three channels IJ=00,20 and 11 are also fitted 
well. Similar situation happens in the IJ=00 channel. But when we try to 
cancel the bound states both in IJ=00 and IJ=20 channels simultaneously, 
it is difficult to find a set of Zj parameters within reasonable ranges to agree 
well with the experimental phase shifts simultaneously in three channels. We 
ascribe this difficulty to the non-perfectness of the Pade approximation itself. 

The existence of the physical sheet resonances violate causality and must 
be a false prediction of Pade approximation. One may argue that such kind 
of false poles locate very far from the physical region and the region where 
ChPT is valid (inside the circle \s\ « lGeV 2 ), and therefore the Pade 
approximation is still acceptable numerically in phenomenology. However 
a careful analysis just reveals the opposite. As can be seen from tables 1 
and 2 the residues of the PSRs are usually very large that in some cases 
it strongly plagues the prediction of the Pade approximation. To see this 
more clearly, taking the IJ=20 channel for example, we recall the following 
dispersion relations 
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for IJ=20, and 



8141 + 1498^ + 2212Z 2 - 1260Z 3 + 454iv / 71n( 



-^±4=) + 10561n 2 (^±^L) = (17) 



for IJ=00. 
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These two equations are obtained by assuming no pole exists on both sheets. 
But pole contributions from different sheets can be added to the above equa- 
tions easily. The left hand integrals can be evaluated by using the expressions 
obtained from the Pade approximant Tq . The left hand side of Eqs. (ITJf) and 
(fl9l ) can be directly evaluated from the Pade amplitudes in the physical re- 
gion. The right hand side of Eqs. ( |T8D and ( ]T9| ) separate different type of 
contributions, i.e., from cuts, resonances, bound states, virtual states or phys- 
ical sheet resonances. It is clearly seen in Fig. that the contribution from 
the left hand integral in Eq. (|18|) deviates far from the experimental value. 
Only after adding the physical sheet resonance in Eq. (pD may we reproduce 
the phase in the physical region. Though in Eq. (|19|) the left hand integral 
works rather well numerically, the predictions of Pade approximants on poles 
and cuts become no longer trustworthy quantitatively, in the situation when 
there exist physical sheet resonances with large residues. 

From the above discussions we realize that the existence of the physical 
sheet resonance plagues the predictive power of the Pade approximant Tq 2 , 
at quantitative level. Similar things happen also in the IJ=00 channel. f] 
Fortunately Pade approximants still work qualitatively in the sense that it 
correctly predicts the existence of the physical resonances, like a and p. In 
here it is worth emphasizing that the virtual state VS2 in the IJ=20 channel 
may be physical. The residue of such pole, though very small, is much larger 



than the residues of BSi and VSj, and it improves the r.h.s. of Eq. (|T8|). 
There is a simple reason for the existence of such a pole which does not 
rely on the detailed form of the Pade approximant Tq 2 . To understand this, 
we recall that around the bound state pole at s = sb the S matrix can 
be expanded as r/(s — s B ) + C + 0((s — sb)) where C is a constant. The 
accompanying zero will occur at sq ~ s B — r/C for any non zero value of C. 
This simple mechanism explains the pair production of BSi and VSi. The 
reason for the generation of VS2 is similar. Remember that S = l + 2ipT and 
ip is purely real and negative when < s < Arri^. When s — > from positive 
side ip — ► —00 therefore if T(0 + ) is positive then 5(0+) — > —00. But since 
S(4m^) = 1 therefore the S matrix must go through a zero (at least once) 
located somewhere between s = 0+ and s = 4m^. It happens that the [1,1] 
Pade approximation gives a negative value to T^ 1 ' 1 \o + ) and T 1 1 ^ 1,1 '(0 + ), but a 
positive value to Tq' 1 ' 1 ' (0+). From this analysis we realize that the existence of 
the virtual state VS2 is a consequence of more general situation (T(0+) > 0) 

6 In the IJ=00 channel the PSR contribution to sin^Tr) becomes also sizable. 
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than that of the concrete form of the Pade approximant Tq' 1,1 ^. In here the 
signs of the Pade amplitudes Tq^~' 1 \o + ), T^ 1 ' 1 \o + ) and Tq^'^(0 + ) coincide 
with the signs of the lowest order amplitudes Tq 2 (0 + ), T^ 2 (0 + ) and T 2 2 (0 + ), 
respectively, and the lowest order amplitudes are unambiguous predictions 
from current algebra. The coincidence of the sign simply reflects the fact 
that the lowest order amplitudes dominate at s ~ 0, as can be clearly seen 
in Fig. H Though the virtual state may really exist, its effect may be very 
small due to its small coupling, as indicated by tables 1 and 2. 

To conclude, the dispersion relations set up in Refs. [I1J and |18| en- 



able us to examine explicitly contributions from different types of dynamical 
singularities - the resonances, the left-hand cuts, and the bound states or 
virtual states - to the phase shifts. Hence a critical examination on the Pade 
approximation becomes possible. We find that even though the Pade approx- 
imation can give a reasonable global fit to the tttt scattering phase shifts, the 
contribution to the phase shifts can be largely from disastrous physical sheet 
resonances in some situations. In such cases the other predictions of the Pade 
amplitudes become unreliable either, at least at quantitative level. However 
we argue that in the IJ=20 channel the virtual state close to s = as pre- 
dicted by the Pade approximation is a consequence of more general conditions 
and seems to exist physically. 
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Figure 2: Various contributions to cos(2<5^ ) and sin(25^ ) in the IJ=20 
channel. Solid lines represent only the left hand integral contributions in 
Eq. fll8|) and Eq. (flQf) . In the dashed line the virtual state (VS2) contribution 
is added; in the dotted line the contributions from both the virtual state 
and the physical sheet resonance are added. The dotted line numerically 
coincides with the dot-dashed line as it should be, the latter is obtained by 
evaluating the Pade amplitude in the physical region, using the k parameters 
from Ref. §. 
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Figure 3: T 2 2 and T 2 4 amplitudes using l { parameters obtained from our 
global fit to the Pade amplitudes in IJ=00, 20 and 11 channels. The bound 
state pole locates at the point where the two amplitudes meet. 
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